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ABSTRACT 


The State Rate Feedback Implicit Model Follower control 
concept is examined in detail from a classical and modern 
control theory viewpoint. State Rate Feedback Implicit Model 
Following (SRFIMF) is a concept whereby control of the dynamic 
resoonse of a system is achieved by the measurement and feed- 
back of a state rate, normally acceleration. In addition to 
a basic description of the concept, emohasis is vlaced on the 
effect of noise in the measurement of the required feedback 
quantities. Control of the pitch attitude of the AV-8A Harrier 
VTOL aircraft is used as an example of the application of the 
control concepvt. The model of the Harrier used in this study 
includes the effect of both sensor measurement errors and gust 


load inputs. 
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ie ROD iON 


A State Rate Feedback Implicit Model Following (SRFIMF) 
Flight controller has been proposed as a possible approach to 
improving the handling qualities of Vertical and Short Field 
Take Off and Landing (VSTOL) aircraft by Merrick at the NASA 
Ames Aeronautical Research Laboratory [l, 2]. The SRFIMF 
concept has potential applications in various types of control 
problems encountered in aircraft design, among which are 
attitude, guidance, and engine control. The concept has not 
been used in actual flight tests, however several detailed 
Simulations [2, 3] have been conducted at NASA Ames ie which 
SRFIMF control was applied. The most notable features of the 
control scheme are that: 

1. The Input-output relationshiv of a system using SRFIMF 
control is insensitive to changes in airframe and propulsion 
dynamic characteristics. 

2. The dynamic relationship of the output to the input is 
approximately that of a second order system whose frequency and 
damping is chosen by the designer. 

3. The system is self trimming and the commanded output 
variable is independent of external disturbances. 

4. The system has good gust alleviation. 

This study presents a detailed analysis of the SRFIMF control 


concept as applied to the attitude control of VSTOL aircraft. 


Hb 





The following discussion is given to clarify the meaning 
of State Rate Feedback Implicit Model Following as used in 
this study. Model following refers to the ability of a con- 
trol scheme to impart specified dynamic characteristics, 
given by a model, to the closed loop system. The model being 
considered here is a second order response in which the param- 
eters of natural frequency and damping ratio are chosen by the 
designer. Typically a second order response is mathematically 
defined in terms of the states position and velocity. The 
state rate of the second order model is the acceleration. The 
State Rate Feedback Implicit Model Following controller, 
studied here, achieves model following by measurement and 
feedback of the system's state rate, acceleration. The result 
is that a priori knowledge of the plant is not required to 
produce model following. 

To illustrate the use of state rate feedback, consider a 
plant of second order. With the states of the system Xy and 
Koy defined as position and velocity, the representation of 


the plant in matrix notation is 


= ote Ute) 1-1 


The acceleration, Kos equation 1s thus 


X. = -cX. -bX. + U(t) 12 


2 1 2 


16 





The values of b and c define the dynamic behavior of the plant. 
Design of a control system would in general require that b and 


c be known. The feedback scheme or control law of the SRFIMF 


ee 


controller is formulated so that the plant dependent quantities 
of b and c are only involved in the total quantity of (~cX,-bX.). 
With this arrangement the quantity, (-cX,-bX,), can be obtained 
by measurement of the state rate, acceleration, minus the 
input, U. Model following by measurement of the state rate, 

as in this illustration, is the basic concept of the SRFIMF 
eeoneroller. 

The intent of this study is to provide a detailed analysis 
of the SRFIMF controller from a modern and classical control 
viewpoint. Particular emphasis will be placed on a basic 
description of the control scheme and the effects of measure- 
ment errors on the output of the closed loop system. The 
first two sections deal with a classical analysis of the 
controller as applied to the attitude position control of a 
general VSTOL aircraft. The third section considers the effect 
of measurement errors on the system from a modern control 
viewpoint. From that analysis, the steady state covariance 
of the state variables, as a result of measurement uncertain- 
ties, will be found. Finally, the previously developed analy- 
Sis technique will be applied to an example where SRFIMF is 
Meeamtor pitch control of the Harrier aircraft. The effect 
of sensor errors will be examined and the response of the 


Harrier to gust inputs will be determined. 


sky, 





The following assumptions are made in this study: 


1. The system is linear. Non-linear effects such as 


control saturation and actuator hysteresis are not considered. 


2. The dynamic response of the plant will be represented 
at a single point by linearized, rigid-body, transfer functions. 

3. Measurement uncertainties are represented as exponen- 
tially correlated noise and the effect of bias error is not 
considered. The basic description and development of the 


SRFIMF controller will be considered in section II. 
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Pie eocR le LLON TOL THE SREP CONTROLLER 


yeeelnHi BASIC SCHEME 

It is desirable for a flight controller to impart 
specific dynamic characteristics to the closed loop response 
of an aircraft system. Often the desired response is given 
in terms of a natural frequency and damping ratio. Piloted 
Simulations at NASA Ames [2] have indicated that a desirable 
response from the pilot's stick to the aircraft attitude is 
a second order response whose natural frequency, Wat 1s 
approximately 2 rad/sec and damping ratio, ¢, of 0.75. One 
possible approach to the design of an attitude controller is 
to apply model following techniques. Mathematically the 
model for attitude control dynamics can be represented by a 
transfer function in the frequency domain. If, as an example 
we let 6(s) revresent the aircraft pitch angle and 6(s) the 
elevator control input, then the transfer function of a model 
for pitch attitude control can be written as 


aS). 1 3-1 


6 (s) 5° “le Zeus “ we 





In this example, if the desired response has a natural freq- 
uency of 2 rad/sec and damping ratio of 0.75, and equation 


2-l becomes 





io 





This transfer function will be used to illustrate model follow- 
ing techniques in the development of the SRFIMF concept. The 
desired model response, given by 2-1, can also be expressed 

in state variable notation with states Z, of position and Zo 
velocity. The matrix of coefficients of the model is given 


the symbol L and equation 2-1 can be written as 


7 eens (BU 
( i | 
ae Z + U 2-3 
ae 2tw ab 


) 


Given somewhat arbitrary open loop plant dynamics, the object 
of model following is to produce a closed loop system whose 
dynamic characteristics are given by equation 2-1 or equation 
2-3. 

For a general system 


ee BU 


2-4 
£ = "CX 


EWemenyecesot amplicit model following is to force the output 


of the system to follow the model equation 


Z= LZ + BU 2-5 


That is to say, the output, Y, should approximate Z (Y=Z) so 


that 
ee DEG 


By proper choice of the feedback gain F, as shown in figure 


Scho woutput Ls forced to folllow equation 2-$. The 
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implied model shown in the upper portion of the figure, is 


not actually generated but implied by the behavior of Y. 


MODEL fi 


L,B ae Uliana I see eA JESS Te ue es te ae 


PLANT 


r A,B 
[F] | 
Led 


Figure 2-1 Implicit Model Following scheme 





The formulation of the feedback law required in model 
following has been determined by Taylor [4] using optimal 
control techniques. In addition, Erzberger [5] has defined 
algebraic methods for determining if perfect model following 
can be achieved. This algorithm will be examined later in 
this section. In aircraft applications certain physical 
facts about the system allow a form of implicit model follow- 
ing to be obtained by simple intuitive reasoning. 

Consider the problem of stabilizing the pitch attitude 
of a VSTOL aircraft. The following assumptions can be made 


and will lead to a simolified SRFIMF controller: 


AI 





1. The desired model response is of second order with 
frequency, Wo! of 2 rad/sec and damving ratio, ct, of 0.75 
from equation 2-1. 

2. Measurements of angular position, angular rate and 
angular acceleration are available. 

peetne COntrol U is either a thrust or control deflection 
whose net result is to produce angular acceleration of the 
menicle. 

4, ‘The open loop plant is arbitrary and the transfer 
function for 6, 8(s)/U(s) = G(s), may be unknown. 

5. The control law which produces model following will 
be developed so that it represents the difference between 
the vehicle's angular acceleration, as measured, and the 
angular acceleration which would be implied by the model given 
iw 2, above. 

From assumption 1 we obtain the desired closed loop response 


as that of the model and express it as 





where @(s) is the pitch attitude and the constants K and Ky 


are defined for convenience and will be used throughout this 


study as 
Kam 2 cue 3 2-8 
me n 
K = wy? = 4 2-9 
ea n 


Bie 





By requiring the model, equation 2-7, to hold we observe that 


the implied angular acceleration is 
2 
s"6é(s) = ~K,.8 (s) = K3s6 (s) + U(s) 2-10 


6(s) and s6é(s) are angular position and angular rates which, 
by assumption 2, are available from measurements. The 
quantity U(s) is the input to the open loop plant. Given 
eas), sSO(s) and U(s), the quantity s“6 (s) can be calculated 
from equation 2-10. 

We will define W(s) as the control law for the system. 
W(s) will be taken as the difference between the implied model 


acceleration 5“8 (s) and the actual measured acceleration 


5“@ (s) ~ (measured acceleration) 2-11 


W(s) 
Using equation 2-10 and substituting for s“6 (s) one has 


W(s) = ~K,.8 (s) ~ Kys6 (s) + U(s) - 


(measured acceleration) 


Further the quantities 6(s) and s6é(s) will be determined by 


sensors so that 


W(s) = -K, : (measured position) - K+ “(measured 
rate) + U(s) - (measured 2-13 


angular acceleration) 


a 





The quantity U(s) could be considered nilot or other control 
input to the plant. Since we also defined W(s) as the control 
there 1S some ambiguity in the notation. The symbol U(s) will 
be used consistently for the control feedback quantity as 


indicated in figure 2-2. 





Figure 2-2 SRFIMF position controller develoved 
EyeIneuLelve acoumMenes. 


Figure 2-2 shows schematically the control law and unknown 
plant G(s). It can also be seen that the olant input, U(s) 
is a feedback quantity and the implied acceleration, 570 (s) 
is compared to the measured acceleration at point 2 of the 
figure. The term state rate feedback 1s apnlied to this tyne 
of control because of acceleration feedback. It is this 
measurement which brings to the controller information about 
the plant making a vriori knowledge of the transfer function 
for @(s) unnecessary. This aspect of the SRFIMF controller 


is most clearly seen by apnlication of the Erzberger criterion 


mentioned earlier. 
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The Erzberger criterion for exact implicit model follow- 
ing 1s obtained by analysis of the system in state space 


representation where the system is given by 


a ee a oad 


Me Cx 4M 3 


For implicit model following the output is approximated by 
Y = LY Rois 
Taking the derivative of equation 2-15, one has 


Y = CX 2a 


Substituting from equation 2-14 


Ke 
ll 


C(AX + BU) 
oie 2-18 
v CAX + CBU 


A 
S 
ll 


but from equations 2-15 and 2-16 one has 

Y= LY = LCX oad 
Bemenat one has, On equating 2-18 and 2-19, 

CAX + CBU = LCX 


or 2-20 
G3use=—(bC — CA): X 
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Bouvang £Or tne control, U, from 2-20 we have that 
U = jice) * * (LC - cal) X Z=2 1 


wnere (CB) is the vseudo, or generalized, inverse of (CB). 
Eliminating U from 2-20 and 2-21, the condition for perfect 


model following becomes 


MGENGS) @ 117 (EC) = (cA) |= x = 6 2-22 
The use of the pseudo inverse is based on the property of the 
pseudo inverse that (CB) (CB)~ 1S an orthogonal vrojection 
operator on the range of (CB). It then follows that 1f£ 2-22 
holds for all X, the range of CB must contain the range of 
(LC-CA). This implies that 2-16 is valid which has already 
been assumed to be the case. 

We will now apply equation 2-22 to the pitch attitude 
controller given earlier. In this case @(s)/U(s) is assumed 
to be the transfer function of an arbitrary second order 
system. 


oats) 1 


U(s) 5° + bs +c 





The desired closed loop performance is given by the model as in 
equation 2-7. For this example we initially take the feedback 
quantities to be angular vosition and angular rate. Defining 
the quantities in state variable representation we have the 
plant where 


X, = 9(t) 


Ao = 0(t} 
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De 
tl 


ee oe 
0 1 xy 0 
= - U 2-23 
-c —-b X5 a 
and 
Y= CX 
or 
J. 0 xy 
Y = 2-24 
0 i x 


Calculation of the pseudo inverse depends upon the relative 
rank of C and B in this case, following Noble [6], the pseudo 


inverse of (CB) is given as 


1 ° 


m(coaee(Cey- ates) 4 


a: 


From 2-23 and 2-24 


(CB) = if 


therefore 
(CB)” = [0 1] 


To determine if perfect model following is possible, we 
Em cemeltewlmto 2722 Using 2-23, 2-24, 2-8, 2-9 and the 


above for (CB) and (CB) 


Zeal 





fies) (CB) = 1] 


— Z—25 
0 0 
0 0 
[(£C) = (CA)] = 2-26 
~ (Ky + Cc) ~ (Ky + b) 
The result is that 
mice(eE) =1) | (c)-(caA) ]) = O 2 aoF 


Equation 2-27 shows that perfect model following is possible 
for ali X using position and rate feedback only. The control 


can be determined by equation 2-21 with the result 


x 
we é 1 iz 
Ue = [ ( Ky + Cm Ky <2 joy)e) A 2-2 
= “KYA, “RY Xs = CX, sr DXo eae 


The control law given in equation 2-29 requires that the 
constants c and b be known in order to produce the desired 
model following. In the previous discussion we stated that 
the addition of acceleration feedback provided the needed 
information about the plant. To show this, note that the 
latter two terms of equation 2-29, CX, os DXo, can be inter- 


preted in terms of the acceleration of the system, Kor since 


ee xX eh xX. + Us) Dat 


2 i Z 


one has 


Xo - U(s) = ~CX) 7 ~bX. yee ale 
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Or rewriting equation 2-29 one has 


Uae = “Kk XxX 


C x 1 ee UNS) ew 


We can compare this to equation 2-13 rewritten in the same 


m@atation 


W(s) = phe Os KeX, + U(s) - (measured 2-33 


angular acceleration) 


Note, Xy and x, are the measured angular position and measured 
angular rate. This comparison shows that the measured accel- 
eration supplies the terms needed for perfect model following 
without a need for a knowledge of the plant dynamics. 

The significance of this result is that state rate feed- 
back can be used to provide information about the plant in 
applications where the plant has unknown and changing dynamic 
characteristics. This conclusion was reached based upon the 
assumption that the plant was of second order. We shall now 
consider a more general form of the SRFIMF controller and we 
will show that model following can be achieved by measurement 


of state rate for a higher order plant. 


B. THE GENERAL FORM OF THE SRFIMF CONTROLLER 

The preceding discussion presented an intuitive descrip- 
tion Of the principal operation of the SRFIMF controller and 
the use of state rate feedback. It is the intent of this 


section to develop a general form of state rate feedback and 


Zo 





to show again that model following can be achieved without 
a priori knowledge of the plant by using measurement of the 
State rate. We begin by examining the basic SRFIMF controller 


as developed by Merrick at NASA Ames. The block diagram of 


this controller is shown in figure 2-3. 


sX(s) 





Figure 2-3 The general form of the SRFIMF 
controller. 

Assume that figure 2-3 represents a velocity controller, then 
the feedback quantities X(s) and sX(s) are velocity and accel- 
eration respectively. It can be seen that at position 1 we 
are summing {- acceleration + u(s)). It will be shown that the 
transfer function between the input X(s) and the output X(s) 
Goes not depend upon the plant transfer function G(s) and 
that the closed loop response is that of a second order system 
whose damping and natural frequency are determined by the 


choice of feedback constants Ky and Ky. 
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From figure 2-3 W(s) is given by 
K 
W(s) = U(s) - sX(s) + [X,(s) - X(s)]+|-~S_] 2-34 
and for V(s) we write 
V(s) = U(s) - sX(s) Za 


from the definition of the transfer function G(s) we have 


that 
X(s) = W(s) G(s) 2~36 
Combining 2-34 and 2-36 we have that 


OS) ee) ee 


X(s) =4 U(s) - SX(s) + = Ke G(s) 2-37 
Since U(s) = W(s) as seen from figure 2-3, we have 
X(s) = U(s)* G(s) 2-38 


and equation 2-37 can be rewritten as 


cae 6 (Sees (5) Gis uct c (See (5)) Ky, 


aE ee G(s) 2-39 
x 


and 


(x, (s) : x(s)] Ky 


sX(s)G(s) ='——-~,———= - G(s) 2-40 
x 


G(s) on both sides of equation 2-40 cancels and the result 


becomes 


K 
aaa == x 2-44 
Ss 
CS s + Kys + Ry 


2) 1L 





Equation 2-41 shows that the closed loop response of the 
SRFIMF is identically that of a second order system whose 
natural frequency and damping ratio is completely determined 


by the constant feedback gains which are given by 





The results of this analysis indicate that the SRFIMF 
controller is an excellent candidate for VSTOL aircraft 
applications. The form of the controller shown in figure 2-2 
could be used to control a position such as attitude while 
the general form shown in figure 2-3 might be used to control 
rates or velocity. It will now be necessary to examine the 
requirements necessary to impliment a SRFIMF controller in 
an aircraft. In particular, the exact relationship between 
U(s) and W(s) must be considered. The next section will 
examine practical examples of SRFIMF controllers ina realistic 


aircraft environment. 


Gs PRACTICAL EXAMPLES OF THE SRFIMF CONTROLLER 

Two types of controllers are illustrated in this section, 
a position controller, figure 2-4, and a velocity controller, 
figure 2-5. A first order actuator 1s included in the plant 
of each controller. Also included is a compensator in the 
control feedback loop involving U(s). The system is not 
realizable without the compensator. These controllers will 
be used in the analysis of this report. 
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Figure 2-4 Physically realizable SRFIMF 
Peston Contino ller - 
V(s) z 
Z 
2 l 
U(s) 
1 (Compens ator) 
A 
cSt 
A(s) 
X(s) 
3 
K + 
Si X SZ 
“ StKy & ° 
es) 





Figure 2-5 Physically realizable SRFIMF 
rate controller. 


ae 





The plant is assumed to be driven by a control actuator 
and because of that a phase lag can be expected to exist 
between the commanded control signal, W(t) and the input of 
eae plant. The effect of the actuator will be modeled by 


using the first order transfer function 





where H(s) is the output of the actuator, W(s) is the input 
and t is the actuator time constant assumed to be 0.1 sec 
for all of the examples considered in this work. Recall from 
a previous section that SRFIMF feedback contained a term U(t) 
which represented the control input to the vlant and that W(t) 
was assumed to equal U(t). Because of the fact that U(t) was 
equal to W(t) the output signal of the controller could be 
used to cancel the U(t) term in the acceleration feedback. 
It was shown that 

acceleration = —CcXy “bX, Se (ee) 2-43 


and it was assumed that 


U(t) = W(t) 2-44 
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Mmiporder for equation 2-44 to be valid, the system would need 
to have instantaneous response to any input. Because of the 
waeton Of the control actuator this is not possible and a 
compensator 1s placed in the feedback loop between W(s) and 

the summing junction at 1, as shown in figures 2-4 and 2-5. 

The compensator transfer function is A(s). A(s) is chosen so 
that the output of the compensator is dynamically identical 

to the input of the volant, after the control actuator. When 
the control is represented by a first order lag, as in this 
case, A(s) is the same as the transfer function of the actuator. 
In cases of higher order actuator dynamics the result is more 
complicated. An algorithm for determining the necessary 
transfer function, A(s), which must be used in the feedback 

is given by Merrick [2]. We will assume a first order actuator 
and a compensator of the form 1/(ts + 1) as shown in figures 
2-4 and 2-5, for the remainder of this work. 

Actual controllers used at NASA Ames are illustrated in 
figures 2-6 and 2-7. Figure 2-6 1s a speed controller. The 
quantities VD and VDD represent the measured speed and accel- 
acceleration. VC is the commanded speed and W(s) is the 
output of the controller. Figure 2-7 is a SRFIMF position 
controller used to control the vitch attitude, 60, of the RTA 
vehicle. The limiters seen in the figures are not included 
in the model studied here. The purpose of the limiter is to 
prevent the control feedback loop from acting as an integrater 


when the plan control is saturated. This condition could occur 
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when a large difference between the commanded variable, X, 

and the input variable, Xr existed and the plant control 

was saturated. In this condition W(s) would increase with- 
out bound and a reversal of the input would be delayed because 
of the very high value of W(t) at the time the control reversal 
was applied. The limiter was not considered here because of 


the assumption of linearlity. 
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Figure 2-6 SRFIMF controller employed on the RTA 
for speed control. 


~~ 
» 





| TRIM Ineur 
ste. ; Aha | = 


VID 


Figure 2-7 SRFIMF controller employed on RTA for 
pitch attitude control. 
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Tit. TIME AND FREQUENCY DOMAIN ANALYSIS OF THE SRFIMF 


A. REPRESENTATIVE AIRCRAFT PLANT DYNAMICS 

The State Rate Feedback Implicit Model Follower controller 
presented in section II, figure 2-4 will now be examined using 
Root Locus, Bode Analysis and Time Simulations. Root Locus 
analysis will show the asymptotic behavior of the closed loop 
poles and Bode analysis will show the filter characteristics 
of the controller. Simulated time response of the controller 
will be shown in order to demonstrate the ability of the 
SRFIMF to follow the given model and to show graphically the 
Poeentrimning f£eature of the controller. The model and repre- 
sentative plant will next be defined. 

The model for the system is taken as a second order system 


with a transfer function given by 


M(s) x 





The gains Ky and K are selected to yield two different models, 
one a position controller and the other a rate controller as 


shown in the following table. 








Pole 
; ane 5 Me = Location 
Betitude Control 2 rad/sec +75 4 3 -1.5 leon 
Rate Control 1.23 rad/ 77 O23 ok 176 -.88 ~891 


TABLE 3-1 Assumed Model Constants. 
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The values listed in table 3-1 are based upon piloted simula- 
tions at NASA Ames [2] and reflect what has been judged to be 
Representative of good handling qualities. 

The plant transfer function is also taken as a second 


order system 


G(s) = a 3-2 


Se + bse 
While initially this assumption may seem unrepresentative of 
mmpeatyatrcraft, proper choice of b and ¢ can represent the 
dynamic behavior of the individual modes of any aircraft, 
provided that first order modes are taken two at a time. To 


illustrate, consider the general transfer function 


M 
K IT (s + Z.) 
1 ae 


C(s) = To = Se ee 3S8 


(s + Ps) I (s¢ + 2t_ wis + ww) 
K=l K K K 


Y 
ay 
ela 


If the poles are distinct then equation 3-3 can be expanded 


into partial fractions as follows 


g is Z 
a b.(s + TZ.,w.,) +C_.W 1-z 
C(s) =5+ - --; is 2 5 ~— —- - 
J=] J K=l Ss + 25, WES a. Wy 
3-4 


It can be seen from the above that the response of a higher 
order system is composed of a summation of first and second 


order terms. When the first order terms are taken two at a 
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time, the total system resvonse can be expressed as the sum 
of only second order terms with constant numerators. The 
assumption of linearity allows us to analyze the modes 
independently. It remains to choose values of b and c which 
represent typical aircraft modes. To do this, it will be 
necessary to survey several representative aircraft. From 
this survey a number of modes will be selected as typical of 
VSTOL applications. It will be assumed that acceptable control 
of each of the modes will lead to accenvtable control of the 
system. This method of analysis leads to considerable 
Simplification since it allows the plant to be taken as a 
second order system. While this simplified plant will serve 
the majority of the analysis, a more complex plant will be 
examined in section V. 

Table 3-2 is a listing of pole locations for selected 
aircraft. Figure 3-1 is a sketch of these poles in the complex 
plane and an assumed envelope of VSTOL pole locations. Table 
3-3 lists points which will be considered by time simulation. 
From this list, four representative modes will be examined by 


moot Locus. 
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AIRCRAFT 


F106 


Lateral/Directional 


Aa 


Longitudinal 


AS 
Lateral/Directional 
VZ4 


Longitudinal 


VZ4 
Lateral/Directional 


H19 
(helo) 


Harrier 
Longitudinal 
Harrier 


Lateral/Directional 


Table 3-2 Typical airc 
characteristic equa 


2 
(s + 1/T,) (s + 1/T.,) (s + 20w Ss + ow, 


FLIGHT CONDITION 


S.L. M¥.2 
20,000', M?.9 


S.L. M2.2 
short period 


long period 


15,000', M=.9 
short period 


long period 
S.L. M=.2 
15,000", Me.9 
QO kts. 

26.5 kts. 


75.00% tS. 
short period 


Tong period 

Oo kts. 

Sek tS. 

0 kts. 

115 kts. 

0 kts. 

60 kts. 

0 kts. 

60 kts. 

120 kts. 


a ee eee 
wy tC 
2.42 62 
3.01  .159 
1.56 .31 
1.52 .087 
623.344 
aess073 
1.89  .05 
6.61  .096 
7a 809 
onicmer oan 
3.4 .374 
316.346 
669 -.347 
1.59 421 
fas e250 
38 0 =.043 
ieee ae 
a2 ot 
150) s50 
oes 
1@  =.15 


Re, 


-1.5 
-.48 


-.48 


-.013 
=.20 


+01 
-.09 
-.641 
tle 
- 085 
-.127 


=~ Lt 
+, 30 
-.67 
+.10 
+.016 
+,148 
+.30 
+.26 
+2350 
+. cr 


OSCILLATORY POLES 


Im. 


eo 
3.0 


1.5 


15 
so85 


12 
1.89 
6.6 
66 
20 
3.4 


. 30 
103 
1.4 
42 
40 
es 
ae 
45 
le) 
1.78 


Source: References [4 & 6] 


Al 


REAL POLES 
1/T, V/T, 
-. 169 -.59 
-1.84 +.006 
-.065 - .56 
-2.48 - .006 
-.82 =F hay 
-.65 -.90 
=.09 -.87 
-.9 -1.05 
-.33 ~.02 
-.073 -1.0 
-.015 -.58 
-.068 -1.26 
-.056 -1.73 


raft pole locations where the 
tion has the form 


2, = 0 





OHW | 


S° 


4 





“SUOTREDOT oTod 4zerAoATe TeotdAL T-¢€ eanhty 
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Mode i _Real _In. b 
1 ane) 3 a2 
2 4 ius Zee Ore 
3 rer ee OS 
4 a Ola. O38 On2 
> =eeey Ol Ore 
6 02.5 is: -0.6 
7 053 C20 a1 5 
8 Ge 0.14 -0.2 
3 00 020 Cr 
10 = oy, 6 

0.6 


Table 3-3 Aircraft modes used in simulations. 
The modes were choosen to be typical of 
VSTOL aircraft as seen in figure 3-1. 
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B. FREQUENCY RESPONSE AND BODE ANALYSIS 

The object of Bode Analysis is to examine the frequency 
response characteristics of the SRFIMF controller. It will 
be necessary to develop the transfer function of the position 
controller shown in figure 2-4. The plant transfer function 


is given in general as 


eS) CPG 
peste Ws) — ao + bs + C 





As in earlier sections the control actuator transfer function, 


H(s), is assumed to be of the form 


z= a 
SS as au) S58 


The transfer function of the plant is given by G(s), the 
output position, defined as X(s) and the commanded input is 
X,(s). The compensator shown in figure 2-4, and described 


in section 2-3 is 


ae 
(ts + 1) 


A(s) 
From figure 2-4 we can write at point 1, the equation 


V(s) = A(s)W(s) - s*X(s) 


where 5°X (s) is now the acceleration feedback. We write the 


control guantity, W(s), from figure 2-4 as 


W(s) = U(s) - K¢sX(s) = Ky [-X(s) - Xo (s) J 3-7 
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and also note 

x(s) = K-G(s)-H(s)-W(s) 3~8 
Combining and rearranging the above expressions we have 
that 


Ga Mei mn 


AAS) K, = R Geen 


Z 
+ ° SL 
oa Kes + K,. X(s) 


and further rearrangement leads to the transfer function for 


the innut-output relation 












+ S + K.S + K 3-9 
x x 


It will be convenient to separate the constant plant gain, 
CPG, from the ovlant transfer function. The remaining transfer 
function, defined as G'(s), has a unity constant gain. Stated 


another way 
G(s) = CPG-G'(s) 3-10 
The parameter KRL is thus defined as 
KRL = K-CPG/t 


empstituting for A(s), H({s), and KRL and, rearranging 3-9 we 
are left with 
KRL°K 


eCSe ee T 5 3-12 
x (s) 7 s(arcey) * KRL (s a Kys + Bee) 





Equation 3-12 shows the effect of the SRFIMF gain parameter, 


2 
KRL. As KRL increases, the model term (s” + Ks + Ky), 
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becomes more dominant in the closed loop transfer function. 
Using equation 3-12, we can examine the frequency response of 
the controller. 

Assuming the plant transfer function to be a second order 
system given by G(s), equation 3-2, equation 3-12 can be 


rewritten as 


KRL*K 
x 


~s 
n 





Q 
Mn 
at 


(KRL + b)s* + (KRL Ke + c)s + K,-KRL 3-13 


The dominant effect of KRL in equation 3-13 is again seen in 
the coefficients of the s and s* terms. From table 3-2 we see 
that a likely range in the values of b and c are -0.6sbs<1.6 
ana 0.0ScS8.64. If KRL is of the order of 25 and b and c are 
in the range of 2 and 8 resvectively, then KRL will dominate 
the terms of equation 3-13 which contain b and c. Equation 
3-13 can be simplified by neglecting b and c with resulting 


transfer function for the closed loop system written as 





X(s) _ RIO. 
= 
XS) 5° + KRL(s~ + KRL Kys “F ) 3-14 


From this later representation it is clear that the plant 
dynamics which were determined by the coefficients b and c no 
longer olay a role in the frequency response analysis. 


Equation 3-14 exvressed in Bode form becomes 


KRL’ RY 
So 








3 


o. (iw ) + KRL (iw) 7 + KRL Ky (iw) + KRL K 
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Equation 3-15 was plotted for frequencies from .01 to 100 
and for KRL values of 1, 10, 25, 50 and 100. Figures G2 co 
3-6 are the plots of the frequency response of the closed loop 
system. It can be seen that the frequency response of the 
system is that of a low pass filter with a break frequency of 
around 3.0 and a 40 db/Dec roll off. It is also seen that the 
response is very nearly in phase with the input for frequencies 
lower than 3.0. The results of the frequency response analysis 
show that the controller possesses good frequency response 
‘characteristics from the input to the output in that it has a 
flat response for all input frequencies of interest and little, 
if any phase shift. It has also been shown by Merrick [2] that 
the controller attenuates plant disturbances in the form of 
applied accelerations. The question of control disturbances 
will be examined from a different point of view in section V. 
We shall now consider the requirements for the gain, KRL, by 


root locus analysis. 


ee ROOT LOCUS ANALYSIS 

The intent of the root locus analysis is to determine the 
magnitude of the SRFIMF controller gain parameter, KRL, neces- 
sary for acceptable model following. Plots of the root locus 
of the oscillatory pole will show that the desired closed loop 
pole location is approached asymptotically as KRL is increased. 
Acceptable performance is determined by the specific applica- 
meen Lon which the SRFIMF controller is used. We begin by 
examining the equations for the closed loop system developed 


in the preceeding section. 
47 
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The characteristic equation for the system can be obtained 
by setting the denominator of equation 3-13 to zero. The 


resulting equation is 


KR S” + (c + K, KRL)s + KRL-K, = 0 aenlG 


At least one root of 3-16 must be real and the factor corres- 
ponding to the real root is defined as (s + £). Assuming that 
the system described by the transfer function, equation 3-13, 
represents a model following system, then the second factor of 
3-16 must be given by the model. In other words, assuming 


perfect model following, equation 3-16 must factor as 


2 x a 
(s +£) (s”~ + Kss + K.) = 0 3-17 


As KRL becomes large, the SRFIMF closed loopv system response 
does approach perfect model following, and it can be stated 


that as KRL increases the condition is approached where 


s? + (b + KRL) s7 + (c + K»KRL) s + KRL* KL ——» 


24; 
(s + £):(s° + Kys + K)) Bae 


In other words, the dynamic behavior of the system's oscilla- 
tory mode is approximately given by the model parameters Ks 
and K,. Expanding the right hand side of 3-18 we have for 


increasing KRL 


5° + (b + KRL) $7 + (c + KsKRL)S + KRL-K, —s 


3 2 3-19 
s~ + (Kk, + f)s + (K + Kyf)s + EK) 


x 
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Comparing coefficients reveals that for increasing KRL 


b + KRL > Ky + £ 3-20a 
+ Ke ; = 

Cc KsKRL > Kyf + io 3-20b 

KRL—~}>ft 3-20c 


The relationships of 3-20 are not equalities, but they indicate 
that while KRL is increasing, the difference between KRL and f 
remains finite. 

Therefore, it can be said that one real vole is on the 
negative real axis and its location is approximated by the 
value of KRL since f is apyvroaching KRL as KRL becomes large 
relative to b or c. To examine the behavior of the oscillatory 
pole as KRL is increased, it will be necessary to rearrange 
equation 3-9 into root locus form. Setting the denominator of 
3-9 equal to zero yields 


1 - A(s)]|. 1 


a ee 2 —4 = 
x Secs) +s + Kys oa K,. = 0 3-21 


Substituting for A(s), H(s) and G(s) and rearranging, the 


general relation for the rool locus is 


2 
ee KRL (s~ + Kss + K.) 


= G'(s) 3-22 


where G'(s) is an arbitrary plant transfer function divided 


by its constant gain, CPG, as was shown by equation 3-10. 
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Equation 3-21 will now be used to determine the root locus 
for the position controller where G'(s) is given by one of four 
of the modes assumed to be typical of VSTOL aircraft as listed 


in table 3-4. The cases considered are 


Mode # from 





Table 3-2 b Cc 
Case 1 -.6 Ze16 
Case 2 6 ol 
Case 3 0 0 
Case 4 10 ec 6 


Table 3-4 Plant Parameters for Root Locus Analysis 


The root locus computer program developed by Melsa and 
Jones [7] was used to evaluate and plot the root locus given 
by equation 3-21. The gain constant, KRL, was varied from 
0 to 100. The resulting root locus trajectory of the oscilla- 
tory poles for the four mode cases are shown in figures 3-7 to 
3-10. From the figures it can be seen that in all cases the 
oscillatory pole approaches the desired value given by the 
model, in this case -1.5, 1.321. Because of scaling, the real 
pole described earlier is not shown in the figures. It can 
also be seen that the pole location of the closed loop system 
1s within 5 vercent of the desired value for KRL of between 
25 and 50. Although a value of KRL equal to 25 is normally 
sufficient, a KRL value of 50 will be used for the remainder 


@®e this work. 
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iH 


KRL = 5.0 
REAL POLE = -2.26\ } 


\ 


ZA ee 


Relea 
fee, = 8.0 
cH REAL POLE = -4.0 
| REAL POLE= 
=1..35 
KRL =20 
| 5 yi, REAL POLE = -17.0 
KRL = 14.0 a 
pe KRY = 99 i KRL= 
- REAL POLE = -96.0 
| A Meine 50 
AS 
KRL = 0 


IMF posi- 
Ficure 3-8 Root locus of the oscillatory pole of a SRF a 
ion controller. The open loop plant has two poles at the origin. 
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14 | 


KRL = 5.0 
REAL POLE =-2.41 


agen 


eer a ee 


2. 






[At xe = 25 
REAL POLE =<22- 





—————— KRL = 99.1 


ee Y REAL POLE = -96.7 





REAL POLE = 0 


Figure 3-9 Root locus of the oscillatory pole of a SRFIMF con- 
troller. The open loop plant has a natural frequency of .95 
rad/sec and damping ratio of 0.475. 


58 









pe Re OO) 
REAL POLE = 
-.22 


Rios Le 
REAL POLE =#97.7 


H.4 


0.0 a 
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Figure 3-10 Root locus of the oscillatory pole of a SRFIMF posi- 
tion controller. The open loop plant has real poles at -1.0 and 
ai 
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The root locus behavior of a more complex plant will be 
shown in section V where the SRFIMF controller will be applied 
to a transfer function of the Harrier aircraft. The system 
will next be simulated for each of the sample VSTOL aircraft 


modes listed in table 3-2. 


De ORFIMF SIMULATIONS 

The simulated response of a plant under the control of 
SRFIMF position controller was done for all ten of the typical 
VSTOL aircraft modes listed in table 3-2. The time histories 
were generated using the CSMP program, reference [6]. The oven 


loop plant for each of the ten systems was assumed to be of 





the form 
ss’ + bs +e 


with CPG = 1 and b and c given in table 3-2. 

The model was as assumed in section II, a second order 
with a natural frequency of 2 rad/sec and damning ratio, of 
0.75. The constants Ke and K, were 3 and 4 respectively as 
determined by equations 2-8 and 2-9. The simulation of the 
model response is shown in figure 3-9. The rise time of the 
modeled resvonse is approximately 1.2 seconds. The model 
response has a slight overshot (2.9 percent) with a peak time 
at 2.4 seconds. Settling to within 1 percent occurs immediately 
after the peak response. For the simulations, KRL was taken to 


be 50 and the actuator time constant, t, was 0.10 seconds. The 
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input to the system and to the open loop plant was a unit step 
function, representing Kae at t = 0. A sample listing of the 
CSMP source code 1S given in Appendix A. 

Figures eel to 3-21 show the time histories of the oven 
and closed loop response of each of the ten simulated systems. 
It can be seen from the figures that the closed loop response 
of each system is indistinguisable from the model response. 
Figures 3-1] to 3-21 graphically show that the response of the 
closed loop system iS advroximated by the model response regard- 
less of the plant being controlled. In addition, the self 
trimming feature of the controller is seen. In all cases the 
steady state value of the output of the closed loop system is 


the same as the input, Kos namely l. 


onl 
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IV. SENSOR NOISE ANALYSIS 


fe oRE IME STOCHASTIC MODEL 

The SRFIMF controller relies upon measurements to produce 
model following. We shall now examine the effect of measure- 
ment uncertainty on the operation of the closed loon system. 
In order to accomplish this we will consider a position con- 


troller with a representative second order plant given by 


G(s) = ay a 


S° + bs +c 


The position controller will be augmented with sensor noise 
sources and the observability and controllability of the 
closed loop system with these noise sources will be analyzed. 
Errors in the state variables will be determined by covariance 
analysis uSing the Lyapunov equation. 

The measured quantities that are considered to be contamin- 
ated by sensor noise are the attitude position, attitude rate 
and attitude acceleration measurements. Two types of sensors 
will be considered. One is of high accuracy and typical of 
good quality inertial navigation system measurements. The 
second is of lower accuracy and might be considered typical of 
Strapdown sensors. Sensor errors are assumed to be of two 
types, bias and high frequency. Bias errors are usually small 
relative to measured quantities and they are of constant value. 


Bias errors have little effect on the dynamic behavior of the 


control system and will not be considered in this study. 
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High frequency errors can be modeled in various ways. It 
is assumed that the sensor noise sources can be modeled by 
exponentially correlated noise. Exvonentially correlated 
noise 1s obtained by passing Gaussian white noise through a 
first order shaping filter. A first order filter is repre- 


sented by the differential equation 
E = -1/T 6 + W, (t) 4-2 


where W, (t) is a scalar white, zero mean, Gaussian noise of 
constant strength, »p. uw is chosen so that the steady state 
value produced by the filter is the square of the sensor error 
standard deviation, ae Following Maybeck [9], the required 


value of »p as the input to the filter 1s given by 
t= 20° /T 4-3 
: 


where T is the correlation time of the exponentially correla- 
ted noise. Figure 4-1 is a schematic representation of a 
typical shaping filter, in Laplace transform notation. é&(s) is 


resulting exponentially correlated noise. 


u 1 E(s) 
Ss + ee 


Figure 4-1 Exponentially correlated noise 
shaping filter 
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Phe standard deviation, o, of the sensor error is a function 

of the sensors. The values of o and T. were obtained from a 
study by Analytical Mechanics Association (AMA), reference [10], 
in which an aircraft using a SRFIMF controller was studied. 
Table 4-1 is a list of the standard deviation estimates which 
were obtained from the AMA study. Also based upon the AMA 
study, we assume a value of a to be 10 seconds for both sensor 


models. 


. Strapdown Inertial 
Attitude position, a, 12 Orel 
Attitude rate, oO. 0.5°/sec 0.5°/sec 
Attitude acceleration, oc. ieee 0.1°/s- 


(TY. = 10 seconds in all cases) 


Table 4-1 Sensor error model parameters 


The position controller shown in figure 2-4 is augmented 
with the sensor error sources and the resulting system is 
shown in figure 4-2. 

To analyze the closed loop controller, we first obtain 
the state space representation of the plant given by equation 
4-1. Taking position and velocity as the state variable of 


the plant we write 


Position = Xy 4-4a 
Velocity = Xy = X, 4-4b 
Acceleration = X5 = ~CXy - bX5 “P CPGX, 4-4c 


fe 





X(t) 





CPG 


s +bst+c x(t) 


Figure 4-2 High frequency error model, 
BOstervon Controller s 


The additional state variables required for analysis of the 


controller as shown in figure 4-2 are 


X. = output ©£ the actuator which is also input to 


the plant 

Ky = output of the actuator compensator 

x. = state variable representing the error in the 
velocity measurement 

Xe = state variable representing the error in 


acceleration measurement 


X. = state variable representing the error in 
nosition measurement 
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From the figure we can write the control law Wis), with 


x = 0 as 
matey K(X, + Xo)-Ky(X5 + X5) —(-cX,-bX5+CPGX,+X,) + X, 
meee tc) he) xen CPGX, + X, - K,X, - 
Xe - K Xo ye 


The remaining first order differential equations, obtained 


from figure 4-1, are given as 


X = -1/tX, + K/tW(t) 4-6a 
me = -1/tX, + 1/1W(t) 4-6b 
ae 2 
Xo = 1/T XX. + 20.,/T, 4-6c 
, 2 
Xe = 1/T X¢ + 26/7. 4-6d 
° ae nS 2 
X = 1/T xX. “- 200/T). 4-6e 


Combining equation 4-4, 4-5 and 4-6 and expressing the system 
in matrix form; we have the closed loop system whose inputs 


are the random sensor noise quantities given by 


vas 
iI 


ae ee 


KS 
il 


Ker 


where o Oo 
e Dp’ 


L 


oo and o_, are the standard deviation of the position 


measurement error of the velocity measurement error and of 
acceleration measurement error respectively and [9 ] is the 


matrix of input white noise powers given by 


_ ae 
Op = #4 = 20,77, 


*Matricies A, B, Q and C are shown on following page. 
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Having the matrix representation of the closed loov controller 
being driven by the sensor noise inputs, we can analyze the 
effect of the sensor error by considering the controllability 
and observability of the system with respect to the noise 
inputs and outputs of position, velocity and acceleration. 

The linear control system computer programs, developed by 
Melsa and Jones [7], were first used with the result that the 
system is both observable and controllable. In order to gain 
a better understanding of the controllability and observability 
of the system, a second approach, following a method suggested 
by Bryson [ll], was used. In the later method we decounvle the 
system of equations by diagonalizing the system given by equa- 
tion 4-7. This produces a system of equations in modal 
coordinates. The diagonalization procedure requires that we 
calculate the eigenvalues and eigenvectors of the A matrix. 
The transformation matrix, T, from the original coordinates 
to the new coordinates, is obtained by normalizing, in a com- 
olex sense, the matrix of eigenvectors. The matrix quantities 
corresponding to the modal coordinates are denoted with a 


superscript, *. The system, in modal coordinates is expressed 


as 
X* = A*X* + B*OQ 4-8a 
Ve aes 4-8b 
and the transformation is as follows 
xe = pty 4-9a 
ae = pla 4-9b 
pe = Tt ts 4-9c 
cet = CT a 
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We wiil now apply the above coordinate transformation to 
examine the observability and controllability of the system 
when the plant is a specific second order system. 

For this example we chose one of the representative VSTOL 
second order plants given by mode number 6 of table 3-2. The 


transfer function for this plant is given by 


et 4-10 


SeeeOS ete 2. 16 
We will assume the following additional parameters 


t+ = .l sec, Ky = 4, K.. = 3, K= 5, T = 10 sec 


with these choices the gain parameter KRL is 
KRL = K CPG/rt = 50 


Because the system is linear with respect to the noise inputs 
we can solve the problem in general assuming a unit value for 


the individual noise standard deviations, or 


wnere Snr Sy and o. is the standard deviation of the position, 
velocity and acceleration measurement error. The original 
system and the resulting diagonalized system was determined by 
a fortran computer code called MODAL, listed in Appendix A. 


The results are given 


Xx AX + BO x= Ak eS 
vo CX Wee 


Matricies A, B, C, A*, B* and C* are shown on the following pages. 
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emt ceinstructive to observe the A matrix of the system. 
In the upper right corner, rows and columns l and 2 represent 
the plant. Row 3 is derived from the controller terms, 
K-W(t)/t. These terms dominate the matrix. 

The diagonal elements of A* are the eigenvalues of the 
Original system. All of the eigenvalues are negative, indicat- 
ing that the system is stable. The system has one coupled 
oscillatory mode, -1.56, 1.34 in the upper right corner of A*. 
This mode corresponds to the model and will be referred to as 
the model mode throughout the remainder of this study. The 
eigenvalue -47.48, row 3 column 3 of A*, is the mode which 
represents the output of the controller. We refer to this as 
the controller mode. It pvhysically represents the real pole 
whose location corresponds to the value of KRL as discussed 
[mene root locus analysis of section III-C. Similarly the 
eigenvalue -10.00, row 4 column 4 of A*, is the compensator 
mode. The other three eigenvalues correspond to the noise 
mptcer time constant, -1/T,. From the transformed system, 
the controllability and observability of each of the previously 
mentioned modes with respect to the noise input can be 
determined. 

faeine B~ matrix the Q, above the first column indicates 
that column 1 is the input vector corresponding to velocity 
error measurement input. The other columns are annotated 
similarly. The first two rows of the B* matrix indicate that 


the model and controller modes are affected by all three 
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measurement errors and that acceleration measurement error is 
least Significant. The zeroes in row 4 of the B* matrix 
indicate that the compensator mode is unaffected by measure- 
ment noise. This indicates that the compensator does not 
contribute to the uncertainty of the system. 

The C* matrix has been annotated in a way similar to B* 
so the xe x, and xe correspond to the velocity, acceleration, 
and position error modes respectively. From the C* matrix it 
1s seen that the observation of position, corresponding to the 
Birst row of the matrix, is affected by all three measurement 
errors but that velocity and acceleration are only slightly 
affected. We conclude from the C* matrix that measurement 
errors will affect the position, which is the quantity being 
controlled, but that the dynamic behavior, velocity and accelera- 
tion are only slightly affected. This result is perhaps due to 
the low frequency error model which results from the choise of 
ce equal to 10 seconds. 

From the modal analysis it is concluded that sensor errors 
can affect the model and controller mode of the closed loop 
system and that the errors can be observed in the measurement 
of the position. The exact relationship between the output 
quantities of position, velocity and acceleration will be 
determined in the next section where the variances of state 
variable Xi Xo, and X3 will be found by solving the Lyapunov 


equation for the systen. 
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B. COVARIANCE ANALYSIS 
The object of the covariance analysis is to determine the 


standard deviations of the error in the states xX and X 


eae 3 
of the system developed in the previous section. The standard 
deviation of Xy and X5 are of interest because they represent 
position and velocity tracking accuracy. X4 1s of interest 
because it is the input to the plant. A large deviation in 
X, due to measurement error would mean unnecessarily high 
amounts of control energy lost because of measurement errors. 
The state variable representation of the system developed in 
section 4-1, where the inputs to the system are noise sources 
representing position, velocity and acceleration measurement 
error, will again be used. 

The problem of determining the covariance of the states of 
a system in the presence of disturbances is typically encount- 
ered in the design of optimal estimators, the Kalman filter. 
In the estimator problem the system is expressed in the follow- 
ing way 

X = AX + BU + GW, (t) 4-11 
b@ere A is the matrix representation of the plant, B the control 
input matrix, W, (t) white noise disturbance to the plant states 
with 
E(w. (t)] = 0, E(W,(t), We(t)] =u 

and G the input matrix of the disturbances. In estimation 


oroblems one assumes that measurements of the system are made 
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which can be expressed in terms of the states of the system 


by the relationship 


2. = HX + Vy 


Here N quantities are measured with an expected error of Vir 
In the developement of the SRFIMF controller, sensor 

noise was defined as a state of the system and the A matrix 

was augmented to include these states. The covariance analysis 

assumes quiescent operation of the controller, therefore, the 

only input to the system is white noise of strength, u, given 

by equation 4-3. The white noise acts only to disturb the 


sensor error states. The governing equation for the covariance 


of the states is given by the Lyapunov differential equation 


P = AP + PA’ + coc? 412 


where P is the matrix of the state variable covariances defined 


as 


and Q is as before, the matrix of white noise input whose power 
is wp. We have chosen to rename the matrix G of equation 4-1l 
to B because the white noise sources are considered to be the 
inputs to the system. 

Again a unit value of standard deviations was assumed and, 
because the system is linear with respect to the noise inputs, 


we can apply later the set of sensors. We will consider the 
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same plant as used in section 4-1. The transfer function 
for the plant is given by 


G(s) CPG 


ee CC A-13 
CRS e465 + 2-16 





At time equal to Zero no uncertainties exist in the system, 
therefore the initial condition on P is zero. 

A computer program, Vary, listed in Appendix A, was 
developed to solve the Lyapunov equation. The program uses 
an International Mathematical and Statistical Library (IMSL) 
subroutine called DVERK. This subroutine solves a system of 
first order differential equations using a Runge-Kutta 
Algorithm. The resulting covariance matrix, P, is printed by 
the program as a function of time, at a number of discrete 
times. The diagonal elements of the P matrix are interpreted 
as the square of the standard deviation or RMS value (0%) of 
each of the state variable uncertainties as a result of the 
noise input. The data obtained from the computation of the P 
matrix is plotted in figures 4-3 to 4-5 and tabulated in 
Appendix B, table B-1l. The figures show the standard deviation 
and X, as functions of time from 


(sigma) of the states X X 


ez 5 
zero to ten seconds. When necessary, the steady state value 
is shown by a "+" on the figure. 

Figures 4-3 to 4-5 can be interpreted in the following 


way. Each individual curve represents the contribution to the 


total state variable error as a result of one measurement 
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source. For example, in figure 4-3 we see that a unit value 
of standard deviationerror in position measurement results 
in a standard deviation error in Xy 6f approximately 0.8 at 
Six seconds, and a steady state error of 0.98°. The total 
expected value of the uncertainty in Xy is the sum of the 
errors resulting from position, velocity and acceleration 


measurements as given by 


Where the notation Cis refers to the ith state variable and 

j refers to a noise source. Thus, Sip refers to the standard 
deviation of xy as a result of position measurement error. 
Table 4-2 summarizes the steady state values of the components 


of the state variable errors. 


x Xx 


0 2 3 
Position measurement error .980 . 360 Zee 
Velocity measurement error 135 2 NiO 
Acceleration measurement error 2245 089 0.4453 


Table 4-2 Steady state errors as a result of 
measurement errors 
The actual measurements are made either by inertial 
navigation sensors or by strapdown sensors whose standard 
deviations are given by table 4-1. We can now apply the 
results obtained from the covariance analysis to obtain an 


estimate of the state variable errors for actual measurement 
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cases. For example, the steady state standard deviation of 
the error in Xy when the measurements are made by an inertial 
system 1S given as 


cn a Ip 1p + qe 1y + Ga°1. 4-15 


where ao! or and q, are the standard deviations of the inertial 


sensor measurements given by table 4-1 as 


Gae= 0.1 


0.5°/sec 


1Q 
< 
il 


Gs 0.1°/sec? 


and Tp 1y and o,, are given in table 4-2. The resulting 
expected error in Xy in steady state is found from equation 
4-15 to be 0.49°. Ina similar way, the steady state errors 
an) Xs Xo and X3 for both inertial measurements and strapdown 


measurements can be evaluated. The results are shown in 


table 4-3. 
oo 29 *3 
Inertial sensors Oe 4 0.18°/sec i seg° seco 
Strapdown sensors 1.6° 0.58°/sec 3.46°/sec” 


Table 4-3 Steady state tracking errors produced 

by strapdown and inertial navigation sensors 
The result of the covariance analysis indicates that the 
expected uncertainity in Xy and X, as a result of measurement 
errors are not significant. For example, in the case of 


: . . ‘ 0 
strapdown sensors, the expected error in position, X,,is iG 
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while measurement errors had standard deviation values of 
i, 0.5 /sec and Ie Jace. 

To interpret the significance of the expected value of 
the error in Xa, the value of X4 will be determined when a 1° 
Zaput 1S applied to the controller. Refer to figure 4-2 and, 
for this example, disregard the effect of the control actuator 
at the input to the plant. Assuming that initially the con- 
troller is in steady operation with feedbacks equal to zero, 


at the instant the 1° input is applied the value of X. is 


determined to be 
X, = X_ -K -K = 20 4-16 


Comparing the value of the standard deviation in X, to the 


3 
value of X3 when a l° unit step input 1s applied to the con- 
troller, we conclude that the value of the standard deviation 
eEEOr in X4 is only 17.3 percent of the value of X3 when a l° 
input is applied to the system. 

From both the modal analysis and the variance analysis it 
can be concluded that high quality acceleration measurements 
are not necessary and that errors in the measurement of position 
and velocity are not amplified. The errors in the states are 
most affected by position measurement error. This result might 
have been expected because the controller is attempting to track 
to the commanded position and errors in the position measurement 


should dominate the uncertainty in the other states. Figures 


4-3 thru’ 4-5 show that the uncertainty in the states reaches a 
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near steady state value ina short time, particularly in the 
case of Xo. The rise time of the errors, seen in the figures, 
should not be interpreted as a lag in the response, since the 


figures represent the expected error in each of the states. 


C. ALTERNATE MEASUREMENT SCHEMES 

Two possible methods of reducing the number of measurements 
required by the SRFIMF controller will be considered. The 
first case 1s to measure acceleration and integrate to obtain 
velocity and position. In the second scheme, position and 
velocity will be measured and acceleration will be estimated 
from knowledge of the plant. The second case will be used to 
obtain an estimate of the added uncertainty caused by the 
acceleration measurement. Figure 4-6 is a schematic representa- 
tion of the eee case. As before the state representation of 


the system is obtained. 





1 
fie] Velocity Estimate Me s+ 
Ke 
x 


Position Est. 


5 
Ea se Acceleration 
X6 ? = Measurement 


Figure 4-6 First alternate measurement scheme. Measured 
acceleration and implied position and velocity. 
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The state variables are defined as 


= position 


= velocity 


= output of the control actuator 


= output of the compensator 


= velocity estimate 


= position estimate 


acceleration error 


The control input, from figure 4-5 is written as 


W ( 


Ss) = CX, + bX 


seu Gr AC me dee 


2 3 ee 


3 7 


ae 


Following a development analogous to that used in section 2-3 


we obtain 


x 


tl 
Mooo oO oO 8 


the state representation of the system to be 


= AX + BQ 


0 
=~ [5 
pe ik a 
Dy 
-b 
0 


0 


0 0 0 
CPG 0 0 
-(CPG-K+1)/t K/t ~K-K,/t 
-CPG/t 0 -KX/Tt 
CPG 0 0 
0 0 Hl 
0 0 0 
Q= | *a 
TY 
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We apply the modal computer program to this system using 


the plant defined by equation 4-10 in order to determine the 


observability and controllability. 


modal coordinates is 


X* = A* X* + B*O 


-47.48 Oe 
De -1.56 
oO -1.34 
a O-0 O40 
GEO ono 
0.0 0.0 
on0 OO 
-1.07 
7 
-.20 
ba. =e 1210 
9989.68 
ceo 
27 2 
0.0 -0.03 
Ca = Ono One 
-0.99 -0.24 
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The resulting system in 


020 Oo 

Ore 0.0 

Oro Onro 

Ore Or0) 

0.0 0.0 
-10.0 O70 

Ore -o.1 | 

X* 

0 -O.01 -0.42 
0 0038 0.04 
ao aOR 2 O20 





From the A* matrix we can see that this system has two zero 
elgenvalues meaning that the modes of the system which corres- 
pond to the open loop integrators are neutrally stable. The 
modal control vector B* indicates that the velocity and 
position estimate modes of this system are strongly affected 
by the noise input. In particular, mode five (row 5 of B*) 
which represents the position estimate mode has a control 
coefficient from the noise source four orders of magnitude 
greater than the model mode (rows 2 and 3 of B*). We expect 
that in this case acceleration sensor noise significantly 
affects the performance of the controller. 

Applying the algorithm used in the previous section to 
this case we obtain the covariance estimate. Figure 4-7 is a 


plot of the standard deviation (sigma) of Xi Xo and X4 assuming 


BY .2 
USB Ox, 
nN 
>< 
+ 32.2 
> 
ct 
6 Ox, 
ois .o . 
Ky 
8.2 — 
4.2 2.4 42 5.8 8.2 {4.2 


TIME 


Figure 4-7 Error in Xyr Xor and X, as a result of measurement 


errors in acceleration. 
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a unit standard deviation noise source as the acceleration 
measurement error. Data used to plot figure 4-7 is also 
tabluated in Appendix B, table B-2. It is seen from the 
Figure that acceleration measurement error causes unbounded 
errors in xy and X. Unbounded errors in any of the states 
are unacceptable in the controller and unless we can improve 
the method of estimating the position and velocity, we will 
be required to measure these quantities. A Luenberger observer 
could be used to estimate position and velocity, however the 
design of the observer would requre use of knowledge of the 
plant. We wish to avoid using detailed knowledge of the plant. 
The second alternate measurement scheme is the case of 

measured position and velocity with acceleration estimated 
from the plant parameters b and c. The estimation is given 
by 

Acceleration = -cX, -bxX, + CPG x 4-17 


This scheme is shown in figure 4-8. 





Acceleration 





Estimate 


Figure 4-8 Second alternate measurement scheme, measured posi- 
tion and velocity. Estimated acceleration. 
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Following the same procedures as before we define the state 


variable as 


Xy = position 

X5 = Vien © GC ilney, 

X = output of the control actuator 
Ky = output of the compensator 

Xe = position measurement error 


Xe = velocity measurement error 
From figure 4-8 we obtain the control law W(t) as 
W(t) = CK. Ec!) xX, + ao +b) X, —aePG X + Xy 
eee ec) X, f (Ke Eby xe 4-18 


Assuming the same second order plant given by equation 4-10 


we obtain the modal transformation of the system. 


X = AX + BO 
aren 
X* = A*X* + B*Q 


= C*x* 


K 
| 


iMatricies A, B, C, A*, B* and C* are shown on the following 


pages. 
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The results of the modal analysis are similar to those 
obtained when acceleration was a measured quantity rather 
than calculated as in this case. The model and controller 
modes are controllable from the noise sources. The position 
and velocity outputs (rows 1 and 2 of the C* matrix) contain 
sensor noise terms. 

The covariance analysis results are shown in figures 4-9 
to 4-11 and listed in Appendix B, table B-3. Here it is seen 
that the standard deviation of Xi Xo and X. is slightly less 
than in the case when acceleration was measured. This is 
because the quantity taken to be acceleration does not contain 
the additional error of actual acceleration measurement. In 
this case the system requires only measurements of position 
ame velocity. 

The second alternate measurement scheme can be used to 
compare a state rate feedback control scheme to a state variable 
feedback controller, from the viewpoint of increased uncertainty 
in the state variable resulting from the additional measurement 
of acceleration. To make the comparison, assume that the 
measurements are made by strapdown type sensors whose measure- 
ment errors are given in table 4-1. Table 4-4 lists the total 
uncertainty of Xo Xo and X. for both measurement schemes. It 
can be seen from table 4-4 that the overall uncertainty of Xs 
xX, and X3 is increased when acceleration is measured, as 


compared to the second alternate measurement scheme when 


acceleration is obtained without measurement. In a practical 
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xy Ko x3 
Second Alternate Measurement 0.75 02267 7s 1.6 ‘/s* 
Scheme, Calculated Acceleration 
Measured Acceleration eo 0.5607 Ss 3.46 °/s* 
Percent increase of expected 1133 117% 116% 


omeer 


Table 4-4 Steady state error in X,1 Xo and X, Eons 


measurements made by strapdown sensors. 


system, uncertainty would exist in the knowledge of constants b 
and c and it could be expected that this would result in addi- 
tional uncertainty in the value of the state variables. It 
might be possible that the uncertainty in the knowledge of the 
plant could negate the advantage gained by estimating accelera- 
tion as in the second alternate measurement scheme. 

In conclusion, it has been shown that sensor noise does 
not adversely affect the SRFIMF position controller. It was 
also shown that acceleration measurement alone is not sufficient 
for acceptable operation of the controller. The analysis tech- 
nique used in this section will next be applied to the analysis 
of a SRFIMF controller when the plant is assumed to be the 


longitudinal axis of the Harrier aircraft. 
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V. APPLICATION OF THE SRFIMF CONCEPT TO THE HARRIER AIRCRAFT 


A. PITCH ATTITUDE CONTROL 

The preceeding discussion dealt with the SRFIMF controller 
with the plant assumed to be a general second order system. 
In this section, the closed loop system will be modeled using 
the transfer function of a VSTOL aircraft. For this purpose 
the AV-8A Harrier, a jet-lift type VTOL aircraft, was chosen. 
The stability derivatives and transfer functions were obtained 
from reference [12]. The SRFIMF controller concept will be 
applied to the pitch attitude control of the aircraft and the 
analysis will be similar to that done earlier when the plant 
was assumed to be that of a second order system. As before, 
we will examine the root locus, time response and the effect 
of measurement errors. The effect of gust inputs to the plant 
will also be considered. 

The transfer function between the pilot's stick and the 
pitch attitude of the Harrier, at 60 kts, 100 ft/sec, as given 


by reference [12] is 


9(s) _ 0.25 (s* + .246s + .00756) 


GM Acces meaneseMmoG7ece 4 .00747 oe 
ee 2mm ge 2b) is ot 026) eee 5-1b 


(aa) (e = 1073) (s- - .1864s + . 1024) 


OW 





8(s) is given in radians and 6, (s) in inches of stick 
displacement. The above transfer function does not include 
the effect of the control actuator which we include, as before, 
in the model of the controller. 

The denominator factors of equation 5-1lb indicate that at 
60 kts the Harrier has an unstable oscillatory mode with a 
natural frequency, of 0.32 and damping ratio of -0.91, and that 
the constant plant gain, CPG, is 0.25. The plant, given by 
equation 5-1, is used to compute the root locus for the system 
by applying equation 3-21, developed for the controller in 
section III and rewritten here for convenience 


2 
KRL(s”~ + Kes + K..) 


ier ree cre 6 (sh = 30 Bel 


where G'(s) is the open loop plant given by 


G(s) 


G'(s) = ape 





For the root locus evaluation, KRL was varied from 1 to 100. 
The value of the actuator time constant, t, was chosen to be 
0.1 sec and Ky and K, were 3 and 4 respectively. 
The root locus of the oscillatory pole is given in 

figure 5-1. The trajectory of this pole is similar to those 
given in section III where the plant was assumed to be a second 
order system. As before, the closed loop system oscillatory 
poles are near those of the model (-1l.5, ioc KRG values 


of about 25 and greater. Because of the scale of the figure, 
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Figure 5-1 Root locus of the oscillatory poles of the longitu- 
dinal axis of the Harrier uSing SRFIMF position control. Note 
that the open loop zeros are nearly canceled at KRL = 1.0. 
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the real poles are not plotted, however their values are 
listed on the figure for various values of KRL. At the final 
value of KRL (96.6), the system's real poles were located at 
meeeoc, —.2079%and -.036. The pole at -93.28 corresponds to 
the controller mode and its location is nearly equal to the 
value of KRL. This result was discussed in section three. 
The two other poles located at -.036 and -.2079 cancel the 
open loop zeroes of the plant located at -.036 and -.21 as 
seen in the open loop transfer function, equation 5-lb. Zero 
cancellation by a pole could be very detrimental to the response 
if the open loop zeroes are located in the right half of the 
complex plane, becaue it would be unreasonable to expect per- 
fect pole-zero cancellation. This is a problem typically 
encountered in model following techniques when the plant has 
zeroes in the right half plane. 

From the root locus, figure 5-1, it can be seen that the 
dynamic behavior of the SRFIMF controller is unchanged by the 
introduction of a more complicated plant. We shall now consider 
the time response of the closed loop system. 

The controller and plant system was simulated using the 
CSMP program discussed in section III. The source code is 
listed in Appendix A. Simulation of the Harrier transfer 
function is illustrated using a signal flow graph shown in 
figure 5-2. The outputs of position, velocity and acceleration 
are shown in the figure as part of the overall transfer function 


Simulation. For the simulation, KRL was chosen to be 50, T was 
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Figure 5-2 Harrier longitudinal signal flow graph transfer 
LUNG ELOn SLM La elem. 

0.1 sec, Ky and Ky were 3 and 4. Plots of the time histories 
of the closed and open loon pvosition, o6(t), as a result of 
the simulation, are given in figure 5-3. The figure indicates 
that the closed loop dynamic response is nearly identical to 
that of the model as shown in figure 3-5. 

To evaluate the effect of the sensor noise, the state 
space representation of the system is required. Four states 
represent the plant. Defining Xy and position and X5 as 


velocity we write the transfer function for the Harrier, given 
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fa equation 5-] in matrix notation as 


X 

1 0 i 
X, 7 0 0 
X 0 0 
Ky SEOOr ly eH 06736 


0 
0 x, 0 
1 0 XK ic 
+ 5 
0 1 X Sein 
3 
+.4495 -. 4096] Ba .1443 
5-2 


The control vector of equation 5-2 is determined by the numera- 


tor of equation 5-1 so that the states Xy and X, are position 


and velocity. The algorithm for determining the required 


control vector is given by Ogata, reference [13]. 


The remaining state variables are defined in the following 


way 


Xe ="OUCDUEOOr sene CONCEOL aeceuacor, anopue co the plant 


X¢ = output of the compensator 


X. = state variable 


Xe = state variable 
error 


Xo = state variable 


Figure 5-4 is the schematic 


From the figure the control 


W(t) = “KL A) = KeX, = 


representing velocity sensor error 


representing acceleration sensor 
representing position sensor error 


representation of the system and 


law is obtained as 


X = CPG Xp + X¢ = Kix. ees I ees 


is 








.25(s°+.246s+.00756) 
§ +.489s°-.449s -+.06736s+.00747 


(Plant) 





Figure 5-4 Schematic representation of the AV-8A Harrier pitch 
axis with the SRFIMF controller. Including sensor noise. 


Combining and expressing the system in matrix form, the closed 


loop system is given by 


X = AX + BQ i 
0 1 0 0 0 0 0 0 0 
0 0 1 0 ace 0 0 0 0 
0 0 0 1 -.06115 Q 0 0 0 
.007 48 -.0675 045 -. 4906 - 1443 0 0 0 0 
n= SKKOLt “KKy/t -K/t 0 -(.25K+1)/r K/t “KKs/t -K/t -KK,/t 
“Kilt AY x -1/t 0 -.25/t 0 “Ky /t “l/r hy 
0 Q 0 0 0 0 ahd i 0 0 
0 0 0 0 0 0 0 -1/T. 0 
0 0 0 0 0 0 0 0 ae 
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0 0 0 5 
20 /T 0 0 
0 0 0 a 
a 0 O 0 Q= 0 20°/T 0 
a ey 
0 0 0 ; 
0 0 20 2, | 
il 0 0 Pen a 
0 1 0 
0 0 1 


Having the system, represented in state space, the technique 
used in the previous section can be applied. The general solu- 
tion will be obtained for a unit value of sensor error standard 
deviation and later applied to the specific sensor suits assumed 
earlier. The constants required to obtain the numerical solu- 
mmom are: K = 20>: + = 0.1: Ky = 4. Ks = 3; KRL = K-CPG/r = 50. 


The resulting system is 


X = AX + BOQ 5-5 
ana GX 

X* = A*X* + B*Q 576 
met 


(Matricies A, B, C, A*, B* and C* are shown on the following pages.) 


The result of the modal transformation is that the noise sources 


have about the same controllability and observability as in the 
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earlier cases. It can be seen from the A* matrix that the 
closed loop system has two additional eigenvalues. These 
correspond to the shifted poles of the system. Variance 
analysis was done in the same way as in the earlier examples 
and the results are shown plotted in figures 5-5 through 5-8 


and tabulated in Appendix B, table B-4. Table 5-1 lists the 











steady state values of the expected error in states Xi Xo 
and Xe 
Ay Xo Xe 
Position measurement error 0.93 O36 35.5 
Velocity measurement error O79 Oe 290 
Acceleration measurement error Ore2s 0.089 0296 


Table 5-1 Steady state values of 11 F5 and Oe 
aS adeegesult Of position, svelocity and 


acceleration measurement error. 


Comparison of table 5-1 with table 4-2 indicates that the 
values of the standard deviation of the xX and X5 errors are 
identical in both cases even though the two plants are very 
different. The standard deviation of the plan input error 


(X, here, X, in section IV) is slightly different in the two 


3 
cases. This should be expected because of the difference in 
the value of control gain K, required to yield a value of KRL 
equal to 50 and the large difference between the two plants. 
From the variance analysis, we obtain an estimate of the 
expected errors in Xi. X5 and Xe which result from the sensor 


suit assumed in section IV. These values are listed in table 
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oy Ko Xs 
Inertial sensors 0.49 0.18/sec Mos sec" 
Strapdown sensors 57 5 0.584/sec 6.26/sec* 


Table 5-2 Standard deviation of steady state tracking 
errors in Xie X. and X. as a result of sensor noise in 


the Harrier aircraft. 


The Significance of the total tracking error in Xe is 
determined by comparing the standard deviation of the error 
to the value of Xe when a 1° input is applied to the system. 
In quiescent operation, a 1° step input results in an instan- 


taneous value of Xe of 


2 
ae 


.0175rad °- 1° K, K = 1.4 rad/sec 


2 


80° /sec 


The value of x. given in table 5-2 for strapdown sensors 
whose expected error is in the order of 1° is not large 
compared to a ie input. We conclude that sensor measurement 
errors do not adversely affect the performance of the SRFIMF 


HMeGreier pitch controller. 


B. PITCH ATTITUDE GUST RESPONSE 

The pitch attitude response of the Harrier from a gust 
input will be considered by applying the same type of analysis 
ised Garlier. For this purpose assume that the gust acts as 
an additional, uncontrolled input to the system. The symbol 
Sg is used to denote this input. The gust input will be 


a2 





modeled as an exponentially correlated noise source in the 
Same way that the sensor noise was modeled. The schematic 
representation of an external exponentially correlated noise 
source was shown in figure 4-1. In this case the output of 
the noise source, £é(s), is the gust input to the system and 
W, (s) is the white noise input of strength Ug as given by 
equation 4-12. The standard deviation and correlation time 
for the gust can be determined by examining the Dryden wind 
model used in simulation by NASA. At altitudes of 500 ft. 
and above the Dryden model assumes that the RMS value of 


atmospheric turbulence, Toy! is given as 


Sg = 0.2 Vibes 5 


The corrleation time, To? is determined by a characteristic 


length, Loe divided by the vehicle speed, V, or 


Ty = L /V 


The Dryden model gives the characteristic length of the 


turbulence as 


POO feet N= L 


mere mM ts the altitude. Choosing 500 ft. as the flight 
altitude, a flight speed of 60 kts and 15 kts as the value 


of the wind, the standard deviation and correlation time is 


re = ott, Sec, Ty = 6.5 sec 
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The gust input is represented schematically by figure 5-8. 
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Figure 5-8 Gust shaping filter. 


The transfer. function between the gust and position, 
velocity and acceleration of the vehicle can be determined 
from the stability derivatives given by reference [14]. For 
This example we chose pitch attitude position control as was 
done in the previous example therefore the denominator of the 
gust transfer function is the same as the denominator of 
equation 5-1. The numerator is given by McRuer [14] for a 


vertical gust as 
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= é z Z 
= (Mo - Ms f (M Zo (Ma Ms + 
x, (M + ZMq) ae (V-M, + Z,Mq) 5-5 


The stability derivatives for the Harrier at 60 kts are (6) 
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From equation 5-5 and 5-1, we obtain the transfer function 


for the pitch attitude from a gust input 


9 (s) me Gems - Heelan noi 


gt) a>) 4896s° | .4495° + 


-063 


6s + .00747 


inewtzansfer function from the pilot*s stick to the pitch 


angle is given by equation 5-1 where the input to the plant 


is 6: Expressing the plant, with multiple inputs of stick and 


gust load, in matrix notation with states X 


position and velocity we have 


xy 0 il 0 0 x) 
Xo 0 0 1 0 Xo 
X. = 0 0 0 1 Xo 


<e 


|-.00747 -.06736 +.4495 -.4896] Ux 


1 


ae 


and X 


0 -.0026 
£25 = 00278 
-.061  .00306 
1443 .00345_| 


The B matrix has been determined as before so that states 


xy and Xo represent position and velocity. 


5 defined as 


[6 


Having the plant transfer function defined by 5-7 we can 


develope the equation to represent the closed loop COne collet. 


The additional state variables required are defined 


Xo = 
Xe = Output of the compensator 
x. = gust input to the plant, 8g 


(sie 


output of the control actuator, input to plant 


Pc Gem: 


ee 


The schematic of the closed loop system is shown in figure 5-9. 
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Figure 5-9 Schematic representation of the AV-8A Pitch 
control using a SRFIMF controller with wind gust input. 


With Xo = 0, the control law W(t) is 


W(t) = ~K LX, Kas ees -25X_, ay Xe + Ky ee OO 6 ct 


x2 


- 00473) X, 


Rearranging in a manner similar to earlier work, we write the 


closed loop equation in matrix notation as 


X = AX + BO 5-8 
0 ] 0 0 0 0 - .0026 
0 0 1 0 25 0 - .00473 
0 0 0 1 -.061 0 .00306 
A = - .00748 -.0675 ~45 -.49 - 1443 0 . 00346 
here ~KK3/t -K/t 0 -(.25*K+l)/t K/t (Ky + .0026+.00473)K/t 
meet Shey -I/r 0 -.25/t Koi (K; + .0026+.00473)/t 
0 0 0 0 0 0 “VT, 
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B = 0 Q= 224/74 | 


The covariance of the states were calculated for the system 

described above as was done in previous examples. The results 

of the covariance analysis was that the RMS values of xy and 

x are very small and quickly reach the steady state values of 
S2.:= 96 = .0017 oc. =0 = .0134 

1 Vs a 

In the case of the gust we are interested in the pitch accelera- 

tion of the vehicle as a result of the gust input. Acceleration 


is given as Xo and from equation 5-8 is 


Aeecetlersatwoms=— 9 <. + .25X%2)— 2 00473x 5-9 


3 5 7 


Recalling that the elements of the covariance matrix, P, are 
the squares of the standard deviations of the state variables, 


we write from equation 5-9 that 


. , 
55 .00473P 


a, aes 
Oo = P2 + .25P 77 


accel 3 
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ghe standard deviation in pitch acceleration as a result of 
a gust input was plotted along with the gust input and the 
RMS value of X- as a function of time in figure 5-7. The 
results are also listed in Appendix B table B-5. The gust 
is shown for comparison of response time. 

It 1s seen from the figure that the rise time in the 
RMS value of the vehicle acceleration is much slower than the 
gust itself. This result can be intervretated as a smoothing 
of the gust by the aircraft. This is primarily due to the 
slow response of the aircraft to the gust. The controller 
has a relatively fast reaction time and can maintain the 
output xy with only very small errors as a result of a gust 
put . 

The analysis of the SRFIMF controller applied to the 
Harrier has indicated that the qualities of the controller 
found by analysis when the plant is assumed to be a second 
order system apply equally well when the plant is an actual 
VSTOL aircraft. We have also seen that the controller does 
not produce undesireable dynamic response when the vehicle 
is subjected to gust inputs and that the dynamic: response of 
the aircraft and controller is due to the response of the 


meerart alone to the gust. 
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Viz CONCLUSIONS 


The use of State Rate Feedback Implicit Model Following, 
for the attitude control of VSTOL aircraft, has been studied. 
The SRFIMF control scheme is simple, easy to impvlement and 
can be expected to be reliable. It was shown that the 
dynamic response of the closed loop system is a second order 
response and that the natural frequency and damvoing are free 
choices of the designer. It was also shown that sensor noise 
does not adversely affect the operation of the system. Detailed 
conclusions of this study are: 

1. Model following was shown in section III-D by simula- 
tions to be very good. Despite changing plant dynamics, 
representing a variety of flight conditions, the closed loop 
system had a response characteristic of that of the model 
and that the output of the closed loop system was driven to 
the value of the input. 

2. The frequency response of the system is that of a low 
pass filter and that there was no phase shift between the 
output and the input at low frequencies. 

3. Non-minimum phase system, pole-zero cancellation can 
lead to unstable performance. 

4. High quality sensors are not required. It was shown 
by the covariance analysis in section IV-8 that the error in 


position is approximately equal to the error in bosieLon 
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measurement olus 75 nercent of the error in velocity measure- 
ment plus 25 percent of the error in acceleration measurement. 
The general result is that sensor errors do not adversely 
affect the performance of the system. 

5. Estimation of position and velocity by integration of 
acceleration was shown to lead to unstable verformance because 
of neutrally stable modes resulting from the integration. It 
was shown in section IV-C that acceleration sensor noise 
disturbs the neutrally stable modes with the result that the 
system diverges. It is then concluded that the quantities of 
position and velocity must be measured in order to have a 
stable system which does not require knowledge of the plant. 

6. The system is capable of compensating for gust inputs. 
The analysis of section V-C showed that the error in position 
outout was small as a result of the gust. Secondly, the 
response of the system to the gust is smoothed by the action 


Seethne uncontrolled plant. 
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